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Abstract 

We speculate on whether a certain p-adic stabiUty phenomenon, observed by David 
Robbins empirically for Dodgson condensation, appears in other nonlinear recurrence 
relations that "unexpectedly" produce integer or nearly- integer sequences. We exhibit 
an example (number friezes) where this phenomenon provably occurs. 

This note may be viewed as an addendum to Robbins's note |^ in this volume. Its 
purpose is to speculate on whether the p-adic stability phenomenon that Robbins observed 
empirically for Dodgson condensation appears in other nonlinear recurrence relations that 
"unexpectedly" produce integer or nearly-integer sequences, and to provide an example where 
this provably occurs. 

In order to carry out this speculation, we'll phrase Robbins's observation in a somewhat 
more general framework. For us, a recurrence relation over a field K will consist of a finite 
partially ordered set S plus, for each s G 5*, a rational function fg = Ps/Qs over K in the 
indeterminate vector {xt)t<s- (The restriction to S finite does not concede any generality for 
our purposes: to consider an infinite recurrence, look instead at all of its finite truncations.) 
We also assume (for simplicity) that the partial order on S is generated by the relation in 
which t is less than s if fs is nonconstant as a function of Xt alone. In this case, s E S is 
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minimal for the partial order if and only if fg is a constant function; we thus use the term 
initial interchangeably with "minimal" . 

Before proceeding further, it will be helpful to set up some more notation. For / = {is)ses 
a tuple of nonnegative integers, we write for Hsgs^s"' function g : S ^ K, we 

write for Yls^s 9^^)^' ■ Write P, = J2i^s,ix^ and Qs = J2i^s,ix^, where Qsj and bs,i are 
zero for all but finitely many /, and Pg and Qs have no common polynomial factor. 

Suppose now that K is equipped with a discrete (nonarchimedean) valuation v, e.g., 
K = Q with the p-adic valuation for some prime p. Suppose also that the Pg and Qs are 
normalized so that v{asj) > and v{hsj) > for all s and /, and so that for each s, 

rnm{rn:m.{v{asj),v{hsj)}} = 0. 

Suppose further that there exists a function g : S K such that g{s) = fs{g) for all s E S; 
note that g is unique if it exists, and the only obstruction to its existence is the vanishing of 
Qs for some s. That is, g is the unique solution of the recurrence, and satisfies 

E/ Kir 

for all s & S. 

Now fix a positive integer A^. We denote by * any element of K with v{*) > N; here 
we intend that two different occurrences of * may refer to two different numbers. With this 
convention, we have the following simplification rules: 



(1 + *)(! + *) = ! + * 
(1 + *)/(! + *) = ! + *. 

We also have c* = * whenever v{c) > 0. 

Define an N -perturbation of the recurrence as any function g' : S ^ K such that for each 
s G 5, 

In case s is initial, this yields g'{s) = g{s){l + *); this is the same as saying that v{g'{s) — 
gis))>vigis)) + N. 

The point of this definition is that, in the case K = Qp, g' is a possible result of computing 
fsig') using p-adic floating point numbers with A^-digit mantissas. Specifically, recall from 
[S] that a "p-adic floating point number with an A^-digit mantissa" consists of a pair (a, e), 
where the "mantissa" a is an invertible element of Z/p"Z and the "exponent" e is any integer. 
This pair is used to represent any p-adic number ap'^ such that a is invertible in Zp and the 
image of a under the natural map from Zp to Z/p"Z is a. Hence two numbers r and s admit 
the same representation if and only if r = s{l+p'^u) for some u G Zp, i.e., if v{s/r — l) > N. 

One can then reimagine p-adic floating point arithmetic as being carried out with actual 
p-adic numbers, except that at any point in an arithmetic operation, a gremlin may come 



2 



along and multiply any value by a factor of the form 1 + *. In this interpretation, g'{s) is 
then allowed to be any result of computing fs{g') in the presence of such gremlins. (Note 
that any "gremlin factor" applied after adding two numbers together can be absorbed into 
the gremlin factors by which each summand is multiplied. Also, the reciprocal of a gremlin 
factor is itself a gremlin factor.) 

Given an A^-perturbation g', define its projected precision loss r^lg') at s G S* as 

rs{g') = max{v{Qt{g'))}; 

t<s 

this generalizes the notion of "condensation error" introduced by Robbins. Note that the 
projected precision loss is determined by the computed denominators rather than the actual 
denominators, which would be the v{Qt{g)); these often but do not always coincide. Note 
also that rs{g') = when s is initial (because the only term in the maximum is v{Qs{g')) = 
f (1) = 0), and that rs{g') > rt{g') whenever t < s, i.e., the bound gets larger (i.e., worse) as 
you go along. 

We say that the recurrence exhibits Robbins stability if for any positive integer N, any 
A^-perturbation g', and any s G S', if Tg^g') < N, then 

v{g'{s) - g{s)) >N- r,{g') + min{0, vig{s))}. 

Robbins's conjecture in 0, made on the basis of copious numerical evidence, then essentially 
(but see next paragraph) amounts to the statement that the recurrence obtained from Dodg- 
son condensation of a matrix of indeterminates (indexed by the connected minors) exhibits 
Robbins stability. (Note that the term mm{0, v{g{s))} drops out in Robbins's case because 
v{g{s)) is always nonnegative; this seems to be warranted by experimental evidence, as we 
note at the very end.) 

It may be more accurate to speak here of "weak Robbins stability", as we are actually 
generalizing a slightly restricted version of Robbins's conjecture. That is because Robbins 
permits the "borderline" case rs{g') = N; indeed, the comment in [S' that "a quantity can 
be accurate to zero places" suggests that this permission was deliberate. However, we are 
not entirely sure whether to believe the conjecture in the borderline case, and our proof of 
Robbins stability in the one nontrivial case we can prove (see the Theorem below) does not 
handle the borderline future clarification of this issue would be welcome. 

It may be helpful to recall (or rather, to attempt to reconstruct) some of Robbins's 
motivation for making his original conjecture. The quantity N—rs{g') measures the extent to 
which we can distinguish the denominators we have divided by so far from zero. To the extent 
that we can make this distinction, we expect that Dodgson condensation should continue to 
work; this expectation is formalized in the inequality defining the stability property. 

However, the assertion that A^ — rs{g') measures our ability to distinguish denominators 
from zeroes is only really meaningful if those denominators are integral. This suggests that 
it may not be wise to expect stability for recurrences in which denominators occur in an 
unsystematic fashion; this caution is borne out by a simple example, which we now give. 

Take S = {0,1,2,3,4,5,6,7}, equipped with the ordering that agrees with the usual 
ordering except that and 1 are not comparable, and consider the recurrence over Q given 
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by 

Xo = 5, Xi = -5, Xn = (n = 2, ...,7). 

Xn-2 

The function g in this case takes the values 

6 11 7 40 365 663 



5,-5, 



5' 25' 15' 33' 77 '140' 
Let V denote the 2-adic valuation; then the function g' taking the values 

6 11 793 4040 20365 17463 



^' ^' r ' or ' 



5' 25' 15 ' 33 ' 8723 ' 1601860 
is an A^-perturbation for iV = 6, because 

^ 11/25- (l-2«) 

^ ^ ^ -6/5 

and g'{n) = fn{g') for n = 5,6, 7. The projected precision loss is 

rjig') = max{v{5),v{-5), t;(-6/5), f (11/25), t;(7/15), t;(-40/33)} = 3, 

and w(663/140) = —2, so Robbins stability would predict that 

f (-17463/1601860 - 663/140) > - rj^g') + min{0, f (663/140)} = 6- 3- 2 = 1. 

However, —17463/1601860 — 663/140 = —2661195/560651 has valuation 0, so the recurrence 
does not exhibit Robbins stability. 

As noted before, it is unclear whether one should expect Robbins stability to be exhib- 
ited by recurrences with "unpredictable" denominators. However, there is a wide class of 
recurrences in which denominators either do not occur, or occur in a limited and systematic 
fashion; these are the recurrences which exhibit the "Laurent phenomenon" , in the parlance 
of Fomin and Zelevinsky P . That paper establishes that a number of interesting recurrences 
(like Dodgson condensation) have the following property: if one views the initial constants 
as distinct indeterminates, the noninitial terms turn out to be polynomials in these indeter- 
minates and their inverses. (See [Sj for an online discussion of such recurrences and related 
topics.) 

Among recurrences admitting the Laurent phenomenon, Dodgson condensation is but 
one example, and it seems (to us, anyway) that the unexpected cancellations that contribute 
to the Laurent phenomenon may in the condensation case must have something to do with 
the unexpectedly strong bound on the precision loss predicted by Robbins stability. We thus 
pose the question: do other Laurent recurrences exhibit Robbins stability? 

One can trivially construct many recurrences exhibiting Robbins stability, by considering 
those for which Qg = 1 for all s, so that no divisions are ever performed in the calculation 
and hence rs{g') = for all s E S. In fact, these recurrences have a much stronger property. 
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Proposition. Suppose Qs = 1 for alls. Then for any N -perturbation g' , v{g'{s)—g{s)) > N 
(and hence v{g'{s)) >0) for all s E S. 

Proof. We proceed by induction on s; for s minimal, the desired inequality is given directly 
by the definition of an A^-perturbation, so we assume that s is nonminimal and that 

g'(t) = g{t) + * for all t < s. 

In particular, v{g'(t)) > for alH < s. 

We now begin a second induction to show that (g')^ = 9^ + * for all tuples / of nonnegative 
integers indexed by the set of t G 5* with t < s; this induction will be on the sum of the 
entries of /. If this sum is zero, then the desired equality is the trivially true 1 = 1 + *. 
Otherwise, given a tuple / for which the claim is known for all tuples of smaller sum, choose 
some t at which I has a nonzero component, and let J be the tuple obtained by decreasing 
this component by 1. Then g^ = g^ g{t) and likewise for g\ [g'Y = g^ + * hy the inner 
induction hypothesis, and g'{t) = g{t) + * hy the outer induction hypothesis. These imply 
that g'{t) and {g'Y have nonnegative valuation, and so 

{g'Y = {g'YgXt) 

= {g' + *){g{t) + *) 

= 9-^9{t) + 9{t)*+g^ *+* 

= g-'git) + * 

This completes the inner induction, so we may conclude that [g'Y = g' + * for all /. 
To complete the outer induction, note that 

g'is) - g{s) = ^(a,,/ + *){g'Y - a^jg^ 
I 

= E * - E asj* = * 
I I 

since v{asj) > by hypothesis. □ 

On the other hand, it seems not so easy to establish that Robbins stability is exhibited by 
any recurrences, even ones exhibiting the Laurent phenomenon, in which nontrivial divisions 
take place. However, we have succeeded in doing so in one case, which we now describe; it is 
a form of a recurrence of Conway and Coxeter |L , which we will refer to here as the "number 
frieze" recurrence. 

Fix a positive integer n, and set 

S = {{a, b) eZxZ:0 <a<n, <b <n-a}, 
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with the partial order given by 

(a', b') < (a, b) <^=^ a' < a and b<b'<b + a — a'. 

Choose Co, ... , c„_i G -ft' of nonnegative valuation, and define a recurrence on 5" by 

fm = 1 (0 < 6 < n) 

= Cb (0 < 6 < n - 1) 

I{a,b) = (2 < a < n, 0<6<n-a); 

then g exists and takes values with nonnegative valuations. Indeed, as noted in 0, this 
is basically a special case of Dodgson condensation: the f(a,b) are connected minors of the 
tridiagonal matrix 



/ Co 


1 
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1 


Cl 


1 •• 













1 


C2 



















Cn-2 


1 









•• 


1 


Cfi— 


1/ 



and while one cannot condense this matrix (as some of the other connected minors vanish), 
one can recover the number frieze recurrence by instead condensing the matrix 

W(|i-il)(|i-il+l)/2 i-^j^ 

where t is an indeterminate, then setting t = in the resulting polynomials. 

Theorem. The number frieze recurrence f[a,b) exhibits Robbins stability. 

Note that the proof will actually yield a stronger result, as in the trivial case {Qs = 1 for 
all s): it effectively shows that as long as the projected precision loss is strictly less than N, 
Robbins stability holds even using fixed point arithmetic (i.e., working modulo p^) instead 
of floating point arithmetic. 

Proof. Let g' be an A^-perturbation. (To simplify notation, we write g{a, b) and g'{a, b) 
instead of g{{a, b)) and g'{{a, b)).) We prove by induction on a that as long as r(^a,b)ig') < N, 
we have v{g'{a, b) — g{a, b)) > N — r(^a,b){g') (and hence v{g'{a, b)) > 0, since g{a, b) is known 
to have nonnegative valuation); this gives precisely the Robbins stability bound. 

Before continuing, we introduce another notational convention. Put r = r(^a,b){g')i and 
write Y = Z to mean v{Y ~ Z) > N — r (so in particular any star is congruent to 0). Note 
that the congruences Y = Z and Y' = Z' imply that Y + Z = Y' + Z' always; if Y, Z, Y', Z' 
have nonnegative valuation, the congruences also imply that YY' = ZZ' . Moreover, if 
Y = Z and y, Z both have valuation 0, then Y~^ = Z~^ . 
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We now return to the induction. For a — 0, 1, the desired inequahty holds by default 

because (a, 6) is initial. For a = 2, the denominator of j\a,b) is a;(o,fe+i), and g'{0,b+ 1) = 
g{0,b + 1) + * = 1 + * has valuation 0, so again the desired inequality follows. For a = 3 and 
< 6 < n — 3, we have 

,(3,6) = ^^2,%(2,5+l)-l 



^'(3, b) 



g{l,b+l) 
{l + *)g'{2,b)g'i2,b+l) -{! + *) 
+ (1,6+1) 



by the induction hypothesis, g'{2, b) = g{2, b) + *, g'{2, 6+1) = g{2, 6+1)+*, and g'{l, 6+1) = 
5(1,6+1) + *, so 

,(2,6),(2 6+l)-l + *_ 
' ^ 5(1,6+1) + * 

Since Q(a',b'){g') — 1 for — 0, 1, and since for a' — 2 we have as above Q(a',b'){g') = 1 + *, 
we have 

(a,b')<{a,b) 
= v{Q(a,b){9')) 

^v{g'{l,b+l)). 

Hence (since r < N hy assumption) we have ^'(1, 6 + 1) < A^, yielding v{g'{l, 6 + 1) + *) = 
^(^'(l, 6 + 1)); in particular, ^(^(l, 6 + 1)) = v{g'{l, 6 + 1)) = r. We can now write 

_ 5(2,6)5(2,6+1)-! + * 

9{l,b+l) + * 
_ ((5(2, 6)5(2, 6 + 1) - 1)/5(1, 6 + 1)) + (*/5(l, 6 + 1)) 

l + */5(l,&+l) 
_ 5(3,6) + (*/5(l, 6+1)) 

l + (*/5(l,6+l)) 

= ^(3,6), 

as desired. 

Suppose now that a > 4, ri^a,b)ig') < N, and the induction hypothesis holds for all pairs 
(a', 6') < (a, 6); in particular, we have v{g'{a', 6')) > whenever (a', 6') < (a, 6). To eliminate 
some indices, put 

A = g{a-4,b + 2), 
B = 5(a- 3,6+ 1), C = 5(a-3,6 + 2), 
D = g{a-2,b), E = g{a - 2,b + 1), F = g{a-2,b + 2) 
G = g{a-l,b), H = g{a-l,b + l), 
I = g{a,b) 
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and likewise with primes; note that A, . . . , I all have nonnegative valuation, as do ^4', . . . , if' 
by the induction hypothesis. We then have 

B'C - 1 + * 



E' 



A' + * 



^, ^ D'E' - 1 + * ^, E'F' - 1 



^ G'H' - 1 + * 

+ * ' 

because 51' is an A^-perturbation and v{g'{a',b')) > for a' < a. (More exphcitly, the 
definition of an A^-perturbation implies that E' = {B'C'{1 + *) — (! + *))/{A'{l + *)) and 
the like, but the product of each lettered quantity with a star is again a star.) We also have 
four analogous equations without the primes and stars. Moreover, if (a', b') < (a, b), we have 
> 'i^{a',b'){9') by the way the projected precision loss is defined, so the induction hypothesis 
implies in particular that g'{a', b') = g{a, b); in particular, we have 

A' = A,...,H' = H, 

and we wish to show that /' = /. 

By the induction hypothesis, we have v{E') > 0. liv{E') = 0, then G' = G,H' = H, E' = 
E imply G'H' - 1 + * = GH - 1 and E' + * = E. Since N > r, the congruence E' = E and 
the assumption v{E') — imply v{E) — 0, and so {E' + *)~^ = E~^. Consequently 

G'H' - 1 + * 



I' 



E' + * 
GH-1 



E 

= I 

as desired. 

Since the case v{E) = is okay, we assume hereafter that v{E') > 0; then v{B'C' — 
1 + *) > 0, and hence v{B'C' - 1) > 0. Since v{B') > 0, v{C') > 0, and = 'f;(l) > 
mm{v{B'C'),v{l - B'C')}, this is only possible if v{B') = v{C') = 0. 

We now compute 

G'H' - 1 + * 



E' + * 

{D'E' - 1 + *){E'F' -l + *)-{B' + *)(C" + *)(! + *) 

{B' + *){C' + *){E' + *) 
D'E' E'F' - D'E' - E'F' + 1 - B'C' + * 

B'C'E' + * 
D'E' E'F' - D'E' - E'F' - A'E' + * 

B'C'E' + * 
D'E'F' -D'-F'-A' + {*/E') 
B'C' + {*/E') ■ 
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As before, we have D'E'F' = DEF, D' = D, F' = F, A' = A, and B'C = BC. Moreover, 
from the definition of the projected precision loss, we have 

= , max {v{Q(^a',b'){9'))} 

(a ,b')<{a,b) 
> v{Q(a,b){9')) 

= v{E'), 

and so = 0. 

Since r < N, the facts that v{B'C') = and B'C = BC together imply that v{BC) = 0; 
then the congruence BC = B'C' + {*/E') implies {B'C' + {*/E'))-^ = {BC)-\ This together 
with the previous mentioned congruences and the equation 

_ DEF -D-F-A 
~ BC 

yields /' = /, as desired. □ 

Note that in this example, the precision bound given by Robbins stability is not always 
sharp if one fixes (a, b) and varies over all A^-perturbations. For instance, for = Q with 
the 3-adic valuation, take 

(co,...,C5) = (1,3^-1,-1,1,-11,22). 

For m and sufficiently large (say m > 5 and A^ > 2m), the projected precision loss 
is m (achieved by g{l,l) = —3™), but experiments suggest that v{g'{5,0) — g{5,0)) > 
N — m + 5 always. It would be interesting to find a more precise version of the projected 
precision loss that detects such "localized disruptions" , specifically by relaxing the restriction 
that the bound can only get worse with each successive term. Such a formulation of the 
stabihty phenomenon may even suggest progress towards Robbins's original conjecture or 
generalizations. 

Although all our examples have been recurrences over Q, with v equal to a p-adic valu- 
ation, we have taken care to make our setup more general. In particular, one could use our 
framework to look at Robbins stability in Q(x), with v the x-adic valuation. This might 
serve as a bridge between the Laurent phenomenon and Robbins stability. 

We conclude by mentioning some further experiments the first author has conducted with 
Punyashloka Biswal. Namely, we have been applying Robbins's testing regimen to other 
recurrences exhibiting the "Laurent phenomenon" of compute pairs of A^-perturbations 
using A^-digit p-adic fioating point arithmetic (generating the undetermined p-adic digits at 
random), and compare their difference to the projected precision loss predicted by Robbins 
stability. (This is somewhat easier than comparing one A^-perturbation to the exact solution.) 
Two families of examples we have considered, which both appear to exhibit Robbins stability, 
are the Somos sequences 

Xq = Xi = ■ ■ ■ = Xk-1 = 1, Xn+k = 
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for k = 4,5, 6, 7, and the sequences 

_ x'^+l + CXn+l + d 
— 

given in Example 5.4]. Notably, the latter example seems to require the correction term 
min{0, that we introduced into the definition of Robbins stability. 
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